o 

(N 



(N 



Comparing gravitational waves from nonprecessing and precessing black hole binaries 

in the corotating frame 
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Previous analytic and numerical calculations suggest that, at each instant, the emission from a 
precessing black hole binary closely resembles the emission from a nonprecessing analog. In this 
paper we quantitatively explore the validity and limitations of that correspondence, extracting the 
radiation from a large collection of 191 generic black hole binary merger simulations both in the 
simulation frame and in a corotating frame that tracks precession. To a first approximation, the 
corotating- frame waveforms resemble nonprecessing analogs, based on similarity over a band-limited 
frequency interval defined using a fiducial detector (here, advanced LIGO) and the source's total 
mass M. By restricting attention to masses M G 100, WOOMq, we insure our comparisons are sensi- 
tive only to our simulated late-time inspiral, merger, and ringdown signals. In this mass region, every 
one of our precessing simulations can be fit by some physically similar member of the IMRPhenomB 
phenomenological waveform family to better than 95%; most fit significantly better. The best-fit 
parameters at low and high mass correspond to natural physical limits: the pre-merger orbit and 
post-merger perturbed black hole. Our results suggest that physically-motivated synthetic signals 
can be derived by viewing radiation from suitable nonprecessing binaries in a suitable nonintertial 
reference frame. While a good first approximation, precessing systems have degrees of freedom (i.e., 
the transverse spins) which a nonprecessing simulation cannot reproduce. We quantify the extent to 
which these missing degrees of freedom limit the utility of synthetic precessing signals for detection 
and parameter estimation. 
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I. INTRODUCTION 

Coalescing comparable-mass black hole binaries are 
among the most likely and useful sources of gravitational 
waves for existing and planned gravitational wave detec- 
tors like LIGO [I], Virgo, [2], the Einstein telescope [3], 
and LISA [4] [5]. For sources in a suitable mass range, 
the signal these detectors receive contains significant fea- 
tures from the late-stage, strong-field dynamics of the 
black hole merger. Only full numerical simulations of 
Einstein's equations can provide first-principle models for 
this epoch, including all dynamics and emission [^l [7J. 
Given the large computational cost per simulation, rel- 
atively few well-determined models have been produced. 
Most models thoroughly explored the physics and wave- 
form from nonprecessing binary systems 8-10 . By con- 
trast, relatively few simulations of binaries with more 
generic spins have been published 1 and even fewer have 
provided their gravitational wave signal [TH - [T7] . 

Generic precessing black-hole binaries produce a rich 
multimodal gravitational wave signal during their inspi- 
ral and merger. To simplify the interpretation of these 
signals, several authors have proposed transforming the 
results of numerical simulations, computed in the simu- 
lation frame, into a corotating frame |14U18If2"l"] . Empir- 
ically, a corotating frame demonstrably simplifies (and 
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1 Though many simulations have been performed in the analysis of 

merger recoil kicks, in many cases their gravitational wave signal 

has not been described in detail 6. 11 13 . 



dramatically reduces in number) the modes needed to 
describe the gravitational wave signal, even during and 
after merger. On the other hand, post-Newtonian ex- 
pressions for the strain also naturally decompose into two 
factors: an "instantaneous" factor describing corotating 
emission, transformed by a rotation; see, e.g., [321 [53] ■ A 
corotating frame therefore simplifies the comparison of 
different simulations to one another; to post-Newtonian 
expressions; and to phenomenological or analytic models. 
Additionally, studies of how precession-induced modula- 
tions impact low-mass [23] and high-mass [13] detection 
strategies also naturally express their results in terms 
of corotating frame modes and trajectories. The non- 
precessing search strategies currently being used miss 
precessing signals, roughly in proportion to two factors: 
(a) how much the signal precesses (geometrically) [25] 
and (b) how much the corotating-frame precessing wave- 
form resembles some nonprecessing signal model in the 
search QJ [HI [23] • 

In this paper we compare the gravitational wave signal 
from the merger of a generic precessing black hole binary, 
expressed in a corotating frame, to nonprecessing merger 
waveforms. In this special frame, a precessing merging 
binary emits radiation that closely resembles the signal 
from a nonprecessing analog, for sufficiently short epochs. 
We determine the nonprecessing configuration that best 
fits the corotating-frame's (2, 2) mode, as a function of 
reference mass. Though tuned to one mode, selected ex- 
amples suggest that the best-fit simulation usually re- 
produces multiple modes for a comparable epoch. This 
correspondence suggests a simple kludge to reproduce 
precessing merger signals, proposed directly or implicitly 
by authors who proposed the construction of a corotat- 



ing frame O ITM2"T] . If an appropriate time-dependent 
rotation is known, then a suitable nonprecessing source 
combined with a time dependent rotation generates an 
approximate synthetic precessing waveform. We assess 
this procedure constructively, comparing the line-of-sight 
waveforms generated by a real precessing system and a 
synthetic analog. This procedure can synthesize good ap- 
proximations to short waveforms from precessing black 
hole binaries, using physically-motivated choices for the 
precession rate and inspiral-merger-ringdown signal. 

In the most directly comparable study, Schmidt et al. 
[2"T] applied a similar procedure to a small sample of pre- 
cessing merger waveforms, derived via a post-Newtonian 
approximation. As expected from the functional form 
of the post-Newtonian inspiral signal [53], they found 
that the (post-Newtonian) corotating-frame signal nearly 
matched 2 emission from nonprecessing binaries with 
nearly-identical physical parameters. In short, this study 
provided concrete examples to suggest both that the in- 
spiral signal could be efficiently represented in the coro- 
tating frame using a nonprecessing signal with nearly 
identical physical parameters; by implication, suitable 
hybrid precessing inspiral-merger-ringdown signals follow 
by adjoining precessing post-Newtonian and numerical 
relativity signals in the corotating frame. By contrast, 
our study assesses the similarity between our simulations' 
corotating-frame signal and well-studied models for the 
inspiral and merger of nonprecessing binaries. Specif- 
ically, we employ many (191) generic merger signals; 
adopt a physically-motivated diagnostic to assess wave- 
form similarity; demonstrate that even the merger phase 
of generic precessing signals resembles nonprecessing bi- 
naries, in a corotating frame; and, critically, identify and 
systematic limits to the accuracy of a corotating-frame 
approximation. 

While a good zeroth approximation, this procedure 
does omit physics. Nonprecessing waveforms simply can- 
not self-consistently reproduce features tied to the sys- 
tem's kinematics: the orbital phase versus time; the ring- 
down mode frequencies, set by the final black hole's mass 
and spin; the ringdown mode amplitudes, which can re- 
flect spin-orbit misalignment; et cetera. As a familiar ex- 
ample, in post-Newtonian calculations, time-dependent 
spin-orbit and spin-spin terms must be included in the 
orbital phase and calculated from suitable spin preces- 
sion equations. In general, we find the corotating-frame 
waveform carries additional information (e.g., about the 
transverse spins) that cannot be encoded into a non- 
precessing waveform. The presence of these extra degrees 
of freedom can explain the observationally-relevant dif- 
ferences between corotating-frame and simulation-frame 
results. Over sufficiently long time- and frequency-scales, 



the differences between the corotating frame and simu- 
lation frame become startlingly apparent. In particular, 
early- and late-time waveforms generally resemble differ- 
ent nonprecessing systems. As a result, while corotating 
waveforms approximately resemble nonprecessing modes, 
they do so only for short periods in time and frequency. 
In practice, however, real gravitational wave data analy- 
ses are also limited to a narrow frequency interval and one 
line of sight. Using observationally- motivated diagnostics 
to characterize differences between signals, we find non- 
precessing systems are a surprisingly effective analog of 
generic precessing sources. 

We provide an executive summary and detailed out- 
line in Section [IT] In Section |HI| we review why and how 
we compare simulations using only their corotating-frame 
(2, 2) modes; explain our notation; introduce the simula- 
tions used; and describe the IMRPhenomB model [TU] we 
use in our studies. In Section IV we present our results 
of fitting the precessing system. We present our syn- 
thetic signal and include a discussion of the physics we 



V) Finally, our conclu- 
For interested readers, 



miss in such a signal in Section 
sions are presented in Section [VI 
we include two appendices. Appendix [A] describes the 
previously-developed tools and notation used in this pa- 
per to compare detected gravitational waves and to ex- 
tract a preferred direction from a precessing binary; and 
Appendix [B] which describes the simulations performed 
and numerical tests we adopted to build quantitative con- 
fidence in our results. 



II. EXECUTIVE SUMMARY AND DETAILED 
OUTLINE 

As shown by Figure [T] the pre- and post-merger 
leading-order emission from each simulation is well-fit 
by some nonprecessing inspiral-merger-ringdown model, 
in the corotating frame. The procedure we employ to 
compare two signal models is described in Section |IH[ 
Specifically, for each of the simulations in our sample, 
described at length in Appendix [B] and enumerated in 
Table [I] we have transformed to a corotating frame us- 
ing a suitable time-dependent rotation R(t), then ex- 
tracted the (angular) (/, m) = (2, 2) mode time series 

[R'^ihAt) = J dQR-^Y^ (n)* of the Weyl scalar, 
as projected onto spin- weight —2 harmonics. For each 
possible mass M e 100 — 1000M©, 3 we compare this 
time series with all possible IMRPhenomB phenomeno- 
logical inspiral-merger-ringdown waveforms ^>4p B (A) for 



2 Schmidt et al. 12 lj adopt a completely different definition of 
"match" , based on a single polarization and a white power spec- 
trum. Their quantitative statements cannot be directly com- 
pared with our own. 



3 As noted in the introduction and quantified in the Appendix, 
the duration of our numerical simulations was used to select a 
mass interval, such that only our numerical simulations would 
significantly influence the comparisons we performed. Physically, 
at these masses only the late stages of inspiral and merger would 
be accessible to the advanced LIGO detector, at plausible signal 
amplitudes. 



nonprecessing binaries 10 . The mass range is chosen 
at the low end to avoid the use of hybrid waveforms in 
this study, and at the high end to ensure the inclusion 
of ringdown. Motivated by data analysis, we compare 
signals with a complex "overlap" , maximized over time 
and polarization [261 EZ] ; see Appendix A] for a brief re- 
view. We adopt a fiducial advanced LIGO noise curve 
(zero detuned high power; see [28 ) and perform inte- 
grals over 5 — 2000Hz. 4 Using differential evolution [2"5] . 
a hill-climbing algorithm, we select the best-fitting non- 
precessing (iMPhenomB) mass, spin, and mass ratio for 
each precessing binary merger simulation and mass. The 
best-fit parameters are denoted by Xpb and the best-fit 
match is denoted P m ax,corot: 
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Figure [T] shows the fraction of simulations with match 
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greater than a specified thresh- 
within a few percent of unity; 
and, therefore, a precessing model and its nonprecessing 
analog would appear very similar to a gravitational wave 
detector network with plausible signal amplitude [27; . 



As described in Section IV the best-fitting nonprecess- 
ing parameters depend on precisely which time and fre- 
quency interval we adopt for comparison. For the data- 
analysis-driven diagnostic used here, the best-fitting pa- 
rameters depend on the mass adopted for the precessing 
binary. Roughly speaking, at low mass the best-fitting 
nonprecessing model has similar mass ratio and "effective 
spin" [TCJ1 I5DH52] . When fitting the IMRPhenomB model 
to our data, our results suggest the best-fitting binary 
has the same value of the IMRPhenomB "effective spin" 
parameter xpb! see Section [TV] For illustration, Figure [2] 
shows how xpb depends on simulation mass, for a one- 
parameter family of equal mass binary mergers. 

Conversely, at high mass the signal produced by the 
simulation and any nonprecessing model is dominated 
by quadrupole radiation from the final black hole. In- 
evitably, the best-fitting nonprecessing models must pre- 
dict comparable final states. 

Despite adopting a phenomenologically-motivated 
comparison, our method identifies nontrivial, physically 
significant relationships between precessing and non- 
precessing signals. As demonstrated in Section [iVj the 
best-fit nonprecessing model to the early time (here, low 
mass) signal has the physically-anticipated parameters 
needed to reproduce a common orbital phase evolution: 
similar mass ratio and "effective spin" . For this reason, 



4 For binary masses M < 300Mq(40Hz//; oi1 ,), our results are not 
strongly sensitive to the lower frequency limit adopted for the in- 
tegrand. We adopt fi ow = 5 Hz, to insure only the detector noise 
power spectrum, and not arbitrary choices for fi ow , determines 
our results at high mass. 
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FIG. 1: Precessing binaries resemble nonprecessing 
binaries: Distribution of best fit complex overlap Pmax.corot 
[Eq. |l|] between the IMRPhenomB and our simulated signals of 
mass M. The solid black curve shows the median-probability 
overlap P ma x at each mass: half of our precessing simula- 
tions have a fitting factor above the black solid line at each 
mass; the dotted black curves show the 90% confidence inter- 
val, estimated from our set of 191 simulations. For compari- 
son, the red solid and dotted curves show the median fitting 
factor and 90% confidence interval estimated from 58 non- 
precessing simulations; these simulations are not included in 
the previous list. All calculations are performed by compar- 
ing the two models' (2, 2) modes, each in a suitable corotat- 
ing frame. Our simulations have significantly different initial 
conditions and hence durations: many shorter simulations do 
not have enough data to reliably estimate the waveform and 
hence P for M < 250Mq. Our longer simulations, including 
the nonprecessing simulations, are relatively well-fit by the 
IMRPhenomB model down to ~ 100M© . 



we anticipate the corotating-frame modes can be natu- 
rally extended to arbitrarily early times via "hybridiza- 
tion" with conventional post-Newtonian or effective-one- 
body models for the early inspiral. Even though we only 
fit the (2, 2) mode, detailed followup of a handful of sim- 
ilar systems suggests the best-fitting parameters repro- 
duce several modes simultaneously. Indeed, the striking 
similarity illustrated by the top panel of Figure [4] and 
first noted in Figure 12 of Schmidt et al. [18] was used to 
motivate detailed analysis of the corotating frame. Un- 
fortunately, we only possess continuously-parameterized 
models for one mode from a generic nonprecessing bi- 
nary; we defer a detailed quantitative analysis of multiple 
modes to a future study. 

As noted in Section [Vl this formal similarity in a non- 
inertial frame is directly relevant to physical observers 
in the inertial frame, each limited to a single fixed line 
of sight that does not corotate with the binary. A 
nonprecessing signal, combined with a suitable rotation, 
will generate a reasonable facsimile of a precessing sig- 
nal, in the observer's frame. As a result, even though 
detailed simulations of precessing binaries are compu- 
tationally expensive, a computationally trivial proce- 
dure can generate plausible precessing inspiral-merger- 
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FIG. 2: Best-fit nonprecessing parameters depend on 
precessing parameters and mass, and can differ from 
the simulation parameters: Top panel: For selected simu- 
lations studied in this work, a plot of the best-fitting effective 
spin xpb versus the simulated binary's mass. The best-fit ef- 
fective spin changes as a function of mass: while the initial 
and final state separately resemble a nonprecessing system, 
no single nonprecessing system fits the whole time-dependent 
corotating-frame mode. Instead, the best-fitting parameters 
interpolate (sometimes discontinuously) between the initial 
and final state. Colors indicate the Sq series (black); the Tq 
series (blue); and the short simulations of the T series (red). 
Bottom panel: For all simulations used in this work, the ratio 
Mfit/Msim between the fitted and simulated mass. As indi- 
cated by the shaded region, the best-fitting mass can differ 
by up to ~ f0% from the simulated mass. 



ringdown signals. To demonstrate this agreement, Fig- 
ure [3] uses the best-fit corotating-frame IMRPhenomB pa- 
rameters Xpb derived in the corotating frame; gener- 
ates a time-domain corotating-frame signal consisting of 
only the IMRPhenomB (2, ±2); transforms to the coro- 
tating frame to create a synthetic multimodal signal 
R(t)ip4 PB (t,h); and calculates the complex match P max 
between the simulation-frame (2, 2) modes of the syn- 
thetic signal and the original precessing NR simulation 



FIG. 3: Synthetic precessing waveforms resemble pre- 
cessing waveforms: Top panel: The distribution of Pmax.sim 
[Eq. |2|], the match between 4>4, NR2 2 an d Rip4 PB , a synthetic 
precessing signal derived by applying a known time-dependent 
rotation to a best-fitting IMRPhenomB waveform. This distri- 
bution is effectively identical to that shown in Figure [l] Bot- 
tom panel: For the Sq(1.5, 0.6,45) simulation, a plot of the 
match P max , sim and P m ax,corot [Eq. Q]. As this example il- 
lustrates, except for very high masses P ma x,corot — Pmax.sim- 
To avoid over-weighting the single worst case, which has 
Pnax,corot — 0.95 for all masses, we have explicitly eliminated 
all five resolutions of the strongly-precessing Sq(4,0.6,270,9) 
and the equivalent Tq(4,0.6,90,9) from this comparison. By 
contrast, all resolutions and iterations are included in Figure 
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to which the corotating-frame mode was fit 5 
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At low mass M < 1700M Q , these two quantities agree: 
fmax,corot ^ -P m ax,sim- As expected, at high mass 
these two quantities increasingly disagree. This disagree- 
ment does not reflect intrinsic dissimilarity between non- 



In our investigations, wc constructed a family of synthetic pre- 
cessing signals, unique up to overall mass scale M. The figures 
provided in the text describe the performance of the best-fitting 
member of that one-parameter family. We obtain almost iden- 
tical results if the procedure described in the text is followed 
verbatim. 
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FIG. 4: Corotating signal resembles nonprecessing 
signal : Comparison of \rR~ 1 ip4 l m | (top panel) and 
arg_R~ i/j4 2 2 (bottom panel) for a nonspmning q = 4 sim- 
ulation (dotted) with the corotating waveform from the 
Sq(4,0.6,270,9) simulation (thick). The colored curves cor- 
respond to the (2, ±2) modes (red, solid); the (2,±1) modes 
(blue); and the (3,3) mode (green). A timeshift has been 
applied to align the two simulations. This comparison shows 
several modes from the same nonprecessing simulation repro- 
duce the corotating frame modes from another |18j . 



precessing and precessing signals. Instead, it reflects our 
strict application of the best-fit parameters derived from 
the corotating frame to a slightly different problem: re- 
producing ipiNR m the simulation frame. At high mass, 
slightly different choices for ip^ PB are needed to optimally 
reproduce the NR signal in each frame. 6 In special cases 
like Figure |4l a single nonprecessing model can approxi- 
mate the entire signal. 

Our calculations demonstrate these inexpensive mod- 
els provide a fast, surprisingly accurate model for generic 
precessing signals. Moreover, though not described in de- 
tail here, the rotation operation R(t) can be easily mod- 



6 Though the Fourier transform of R(t)ip4,{t) is well-approximated 
by R(t(f))ip4,(f) at early times, when the rotation operation 
changes slowly, at late times a stationary-phase approximation is 
not sufficiently accurate for our purposes. Equivalently, the ro- 
tation needed at late times oscillates on a timescale comparable 
to the quasinormal mode frequency spacing. 
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FIG. 5: Precessing systems break reflection symmetry 
in the corotating frame: For the Sq(4,0.6,90,9) simula- 
tion provided in [14] , a plot of the reflection-symmetric (solid) 
and reflection-antisymmetric (dotted) parts of the corotating- 
frame I — 2 modes, for m — ±2 (red) and ±1 (blue); see Eq. 
S. This figure suggests that precessing binaries radiate in 
a way unlike any nonprecessing binary, even one viewed in a 
noninertial frame; that this missing physics has a comparable 
effect to higher harmonics like the (2, ±1) modes (e.g., dotted 
red vs solid blue); and that both reflection symmetry break- 
ing and higher harmonics are required to accurately model 
the merger phase of generic precessing binaries. 



eled and fit, particularly in a frame aligned with the to- 
tal angular momentum. In principle, comparisons like 
our own allow fits to functions Ap^(A) and i?(t|A), hence 
producing a "synthetic waveform" R(\)ip4 PB (ApB)- Our 
analysis suggests these "synthetic waveforms" can be 
used as a simple method with which to prototype data 
analysis strategies for precessing, multimodal sources. 
However, because nonprecessing models omit critical 
physics, particularly in the high-mass regime emphasized 
in this paper, we do not recommend high-precision cali- 
bration of these synthetic waveforms, particularly if that 
calibration is limited to a study of only the leading-order 
(2, 2) mode. 

Even when viewed in a rotating frame, a nonprecess- 
ing binary retains its intrinsic symmetries: spins aligned 
with the orbital plane. Without precessing spins to 
break symmetry and source other multipoles, the emit- 
ted radiation has strong symmetries: reflection symmetry 
through the orbital plane, insuring that at every instant 
the binary radiates equal and opposite left- and right- 
handed radiation, in mirror-symmetric directions. By 
contrast, our previous calculations [14] and long expe- 
rience with black hole superkicks demonstrate that pre- 
cessing black holes radiate asymmetrically, emitting pref- 
erentially right-handed or left-handed radiation at any 
instant. Equivalently, even in the corotating frame, pre- 
cessing binaries break reflection symmetry through the 
orbital plane. In Section [V] we use reflection symmetry 
to quantify the extent to which a nonprecessing model 
omits critical physics. To quantify the magnitude to 



which reflection symmetry is broken, we use conjugation 
symmetry to define CP-odd (&;, m (t) = — (— 1) 6j _ m (t)*) 
and CP-even (a^ m (t) = (— l) l ai- m (t)*) parts of the 
corotating-frame Weyl scalar (i? -1 ^): 



[R 1 1p4]l,m = ai,m +h t , 



(3) 



Globally, this specific symmetry C corresponds to a re- 
flection symmetry through the z — plane: 

c^ = rp 4 (Tt-e,<py (4) 



I :m 
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(-2), 



(5) 



Lin 



For a nonprecessing binary, the CP-odd term bi. m = 
0. Using a concrete precessing example, Figure [5] 
shows the CP-even (solid) and CP-odd (dotted) parts of 
the leading-order corotating-frame angular modes ipii m . 
Comparing the red dotted and solid lines, this figure 
demonstrates that after merger, the CP-odd part of the 
dominant mode is significant. In other words, this figure 
shows that, to accurately describe this binary's emission 
versus time and angle, we must include a CP-odd part - 
a part no nonprecessing signal could ever produce! More 
broadly, because a precessing binary has more physics, its 
signal remains more complicated than one from a non- 
precessing binary, even in the corotating frame. 

Figure [6] illustrates the practical implications of these 
asymmetries. This figure shows the match P max ,corot 
for a one-parameter family of high-symmetry simula- 
tions: binaries with mi = to 2 , a\ = ah{9,<j)) 1 and 
a-i = ah(9, 4> + 7r) as a function of <f>. By construction, 
each member of this one-parameter family has an identi- 
cal projection of the spin along the orbital angular mo- 
mentum axis (xpb) and more generally identical total 
spin Si + S2; evolves without precession of the orbital 
plane; and produce black holes whose final mass and spin 
does not change significantly with <\> [Table III. Moreover, 
these binaries produce similar symmetric radiation a 2i 2 
for all orientations <j). To quantify this similarity, Figure 
[6] shows the value of 
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evaluated between one reference simulation and all oth- 
ers, at one reference mass; the high and nearly con- 
stant match demonstrates 02,2 is effectively independent 
of c/). 7 Nonetheless, the total gravitational wave signal 
[R ipi]2,2 = 0,2.2 + &2,2 emitted by these binaries mea- 
surably changes with c/>. These differences can be easily 



7 Similar results occur at all masses. In the interests of brevity, 
we do not plot all of the functions a2,2(<W, even though their 
manifest similarity is immediately apparent to the eye. 
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FIG. 6: Does it matter if we omit the transverse spin?: 

This figure demonstrates that transverse spins have a signif- 
icant effect on the gravitational wave signal, via reflection- 
symmetry-breaking between the (2, ±2) modes. For the 
special-purpose V(a, 6, <j>) family of precessing simulations de- 
scribed in the text and Appendix [B] points show the match 
fmax,corot versus the angle (f>, for the fixed angle 6 = 34°. For 
comparison, this figure uses f/> to indicate the match P max ,NR 
[Eq. (frb] evaluated between the corotating-frame (2, 2) mode 
for one member of that family [F(0.6, 34, 240)] and all others, 
directly demonstrating how much these signals differ from one 
another. For another comparison, this figure uses a to indi- 
cate the match P max ,NR.a [Eq. Mm between the reflection- 
symmetric part of the corotating-frame (2, 2) (02,2) of that 
same simulation and all other angles <f>. The strong variabil- 
ity in the first set (ip) and lack of variability in the second (a) 
suggests that reflection asymmetry is the dominant source of 
variation between different simulations. This figure illustrates 
strong, generic correlations between the precise magnitude 
and direction of the transverse spins; reflection asymmetry 
between the (2 ± 2) modes in the corotating frame; and the 
degree of similarity between the (2, 2) corotating-frame mode 
and conventional non-spinning approximations. 



assessed simply by noting the best-fitting parameters and 
-Pmax.corot recovered when comparing with IMRPhenomB 
vary with mass; as an example of the latter diagnostic, 
see the points in Figure[6J More directly, these differences 
can be demonstrated by computing the overlap between 
different simulations [R~ 1 4 , a]2,2 for the values <f>,4>' ', via 



-PmaxNR = max 



\([R- 1 ^4 NR }2,2\[R'' l ^'NR. 



2,2 



\R 1 ^ iNR \[R' V4^r]2,2 



(7) 



One such example is provided by the markers labeled tp 
in that figure. Figure m\ suggests the asymmetry (62,2) is 
principally responsible both for a significant fraction of 
the mismatch 1 — P max ,corot and for all the fluctuation 
in the match versus <p. This correlation is generic: all 
other generic binaries preferentially radiate asymmetri- 
cally as the spins precess, above and below the instanta- 
neous orbital plane. To a first approximation, the domi- 
nant, symmetric part 02,2 resembles a nonprecessing sig- 
nal; a comparison between 02,2 and IMRPhenomB approx- 
imately determines the best-fitting parameters. Partic- 



ularly during at late times when asymmetries between 
(2, ±2) become particularly significant, however, the in- 
fluence of asymmetry can significantly diminish (or en- 
hance) the similarity between nonprecessing and pro- 
cessing signals - in our context, change the best-fitting 
IMRPhenomB parameters and match P m ax,corot- Broadly 
speaking, our results suggest larger degrees of asymmetry 
between (2, ±2) correlate with a generally larger but spin- 
orientation-dependent mismatch between the corotating- 
frame (2,2) modes and IMRPhenomB. 

To summarize, our study suggests that nonprecessing 
signals, suitably rotated, resemble but cannot reproduce 
with high precision the gravitational wave signal from 
generic precessing binaries. With few exceptions, 8 syn- 
thetic precessing signals generated from nonprecessing 
binaries will not be adequate for high-precision parame- 
ter estimation, unless augmented by new physics (i.e., 
transverse spins) and multiple harmonics. That said, 
given the robust similarity of nonprecessing and pre- 
cessing merger signals over observationally relevant in- 
tervals, we anticipate synthetic precessing signals will be 
extremely useful as plausible signals. We strongly recom- 
mend that existing gravitational wave search strategics 
for high-mass merger signals test how reliably they can 
recover the complex multimodal from precessing binaries 
with M > 100M Q using "synthetic" multimodal inspiral- 
merger-ringdown signals, generated simply by viewing a 
nonprecessing binary in a precessing frame. 

Finally, motivated by previous studies and our own re- 
sults [e.g., Figure p], we anticipate that detailed parame- 
ter estimation of high-mass (M > 100M Q ) binary merg- 
ers will require detailed modeling of multiple modes of 
generically- precessing binaries. A corotating frame and 
a nearly- nonprecessing CP-even approximation can re- 
duce the number of functions to fit; in this mass region, 
however, higher harmonics cannot be neglected. 
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FIG. 7: Corotating modes evolve nearly in phase prior 
to merger: Phases of corotating modes in an absolute scale 
(top panel) and scaled in proportion to their angular order 
4>i,m/\m\ (bottom panel), for the Tq(2,0.6,90) simulation. To 
a good approximation all evolve in phase prior to merger. 
After merger, all angular harmonics shown continue to evolve 
nearly in phase, except the (2, 1) mode. Colors indicate the 
(1,1) modes (red, with solid as I = 2, dotted as I — 3, and 
dashed as I = 4); the (2, 1) mode (blue); the (4,3) mode 
(dashed purple); and the (3,2) mode (dotted purple). 



III. COMPARING THE COROTATING AND 
SIMULATION FRAME 

A. Why compare nonprecessing and precessing 
binaries in the corotating frame? 

As we will discuss at greater length in Section |Vj with 
more degrees of freedom in their underlying kinematics, 
precessing binary mergers never look exactly like a non- 
precessing merger, even in the corotating frame. In the 
corotating frame, the gravitational waves from precess- 
ing binaries break symmetries and carry more informa- 
tion than the signal from the best-fitting nonprecessing 



For special quantities constrained by symmetry, like the open- 
ing angle of the precession cone, parameter estimation should 
be reliable, independent of whether the corotating-frame model 
works well, due both to separation of timescales and due to the 
way polarization modulations enter into the signal 



analog. Nonetheless, the additional degrees of freedom 
are very difficult to excite in a quasicircular inspiral. In 
surprisingly many scenarios, a precessing binary may be 
well-approximated by the emission from a nonprecess- 
ing binary, plus a slowly-changing orientation. Qualita- 
tively speaking, the corotating-frame mode amplitudes 
R Wii m extracted from a precessing binary have sim- 
ilar characteristics as the corresponding modes seen in 
a nonprecessing binary. In the time domain, all of the 
corotating-frame modes i? _1 ?/>4 have smooth phase evo- 
lution and (roughly) smooth amplitude evolution. 9 This 



The corotating-frame \R~ 1 ip^i m \ are weakly modulated during 
the inspiral, in two ways. First, the (I, ±m) modes can increase 
or decrease at the same time (in phase), due to either to residual 
eccentricity or suitable spin precession (for suitable modes). The 
(I, ±m) modes can increase or decrease out of phase, as precess- 
ing spins source symmetry-breaking multipole moments. For ex- 
ample, during inspiral the leading-order asymmetry occurs from 
current quadrupole radiation, sourced by the transverse spins 



smooth and simple behavior in the amplitude and phase 
has already been noted previously [l"8rf2"l"] . As a first 
approximation, we also find corotating modes resemble 
suitable nonprecessing ones; see Figure [4] as an exam- 
ple. Because of this correspondence, the corotating frame 
waveforms roughly satisfy the same symmetries as non- 
precessing binaries. For example, during the inspiral the 
modes evolve in phase (i.e., ex m~ 1 &rgR~ 1 ip4i m ); see Fig- 
ure [7| More generally, all modes are approximately phase 
conjugate (i.e., argi?" 1 ^^ ~ -arg(-l) / i?^ 1 -04; _ m )- 
By contrast, as has repeatedly been demonstrated in the 
literature, the simulation-frame modes are both signifi- 
cantly different and more complicated [H] ITM2T] . Ex- 
cept for nearly-nonprecessing binaries, the substantial 
differences between the corotating-frame and simulation- 
frame waveforms are much more significant than the rel- 
atively small numerical uncertainties in this mass regime. 

In short, the corotating frame provides a significantly 
different, simpler waveform with many approximate sym- 
metries. Previous experience suggests these symmetries 
facilitate interpretation and quantitative calculations for 
precessing merging binaries [IH [21]. Moreover, even 
though observers cannot co-rotate with the binary, an- 
alytic calculations suggest the corotating frame signal 
has direct observational relevance, based on separation of 
timescales. For example, as we highlighted in our exec- 
utive summary, comparisons between the (unobservable) 
corotating-frame signals are directly comparable to com- 
parisons between (observable) simulation-frame quanti- 
ties, l max.corr 



-Pmax.sim [Eqs. ( fTp| ) ; see Figure [3] 



B. Comparing a corotating mode with a 
phenomenological model using a complex overlap 

Following previous studies [TH [3S1 [37] , we use a com- 
plex overlap to assess the observationally-relevant differ- 
ences between two signals a, b 



(a\b) = 2 f 

J — 






(8) 



where each of a, b represent some value for the outgoing 

Weyl scalar ipi 2 ,2 — ( 2 > 2 I^4) = JK.mT^*' see ' e -S-' 
Eq. (24) in [H] and the detailed review in Appendix [A] 
In this work, we apply this diagnostic to pairs of signals 
in the corotating frame [Eq. (nj)] and in the simulation 
frame [Eq. ph]. 



The IMRPhenomB family as a continuous reference 
model 



Though numerous, our numerical simulations of non- 
precessing and precessing binaries only discretely sample 
the model space. For this reason, rather than directly 
compare pairs of simulations, we compare our simulations 
against IMRPhenomB 10, , a model for h 22 as a function of 
total binary mass M, mass ratio r\ = m\ra 2 j(ra\ +?7i2) 2 , 
and a single effective spin parameter xpb'- 

x ± = L ■ (Si/ml + S 2 /m|)/2 (9a) 

Xpb = L ■ (miSi/mf + m 2 S 2 /ml)/M 

= [(l + 6)(S 1 /ml) + (l-6)S 2 /m 2 2 }/2 

= X+ + 5x- (9b) 

where <5 = (toi — m 2 )/M. This model employs a 
physically-motived piecewise-continuous expression for 
^22 (/)j expressed as an amplitude \h 22 \ and phase a.rgh 22 . 
At low frequencies their expressions reproduce con- 
ventional stationary-phase approximations derived from 
post-Newtonian theory; conversely, at high frequencies, 
their amplitude model \h 22 (f)\ has the Lorentzian form 
expected from a quasinormal mode ringdown. Parame- 
ters of this hybrid, phenomenological model were set by 
comparing to numerical simulations and to an extreme- 
mass-ratio limit. This model has been calibrated against 
numerical simulations; includes the effect of aligned 
spins; and has been made publicly-available via the 
lalsimulation toolkit. As this code provides a model 
for the strain (hps) rather than the Weyl scalar, we ex- 
plicitly convert between representations using the Fourier 
transform 



$4pB(f\M>V,XPB,t mgr ) 



1 



(2^/) ; 



;hp B (f) (10) 



being transported around the orbit. The theory underlying the 
latter scenario is briefly reviewed in the text and Appendix [A] 



Because the inspiral and final black hole depend on 
the component masses and spins in distinctly different 
ways, no one IMRPhenomB model can reproduce a generic 
corotating-frame (2, 2) mode for all time. For this reason, 
we distinguish between the simulated binary's param- 
eters (M S i mi r] s i m , xpb, sim) and the parameters of some 
best-fitting IMRPhenomB model (Mg t ,?7fit,XPB,fit). 

By comparing a model against our simulations rather 
than comparable simulations against one another, we 
introduce systematic error. To quantify the degree of 
similarity between the IMRPhenomB model and our non- 
precessing simulations, we have evaluated P m ax,sim = 
-Pmax.corot between this model and an array of 34 non- 
precessing simulations. As indicated by the red curves 
in Figure IT] for most simulations and a wide range of 
masses, some member of the IMRPhenomB model matches 
each of our simulations to better than ~ 1 — 2%. The 
best-fitting parameters are always close to our simulation 
parameters. 

Since IMRPhenomB was first published, other models 
have been developed that may even more accurately re- 
produce the merger of precessing binaries over a wide 



range of masses and spins [5J [33J . While we have com- 
pared these models against some of our numerical simu- 
lations, in this work we limit our comparisons to the suffi- 
ciently accurate and easy-to-evaluate IMRPhenomB model. 



D. In the corotating frame, nonprecessing and 
precessing binaries resemble one another 

Figure [l] summarizes the results of our comparison be- 
tween IMRPhenomB and each numerical simulation, in the 
corotating frame. In short, the two nearly agree: for ev- 
ery simulation, some nonprecessing IMRPhenomB model 
exists which nearly reproduces that simulation's coro- 
tating (2,2) mode ([-R _1 i/>4]2,2) over an observationally 
relevant interval. 

A closer investigation of Figure [TJ however, suggests 
that systematic differences exist between corotating- 
frame waveforms and nonprecessing signals. In that fig- 
ure, the red curves show that truly nonprecessing simula- 
tions are better fit by IMRPhenomB (red curves) than the 
corotating-frame signal from a precessing binary (black 
curves), with matches of order 1 — 3% better. This level 
of disagreement is large enough to be observationally ac- 
cessible: typically, parameter estimation strategies can 
resolve differences between models when their matches 
differ by 1/p 2 , expected to be of order > 1/10 2 for the 
first detections at signal amplitude p ~ 10. 

As highlighted in the executive summary and as we 
will describe in detail later, these differences between 
the corotating-frame and simulation-frame signal are ex- 
pected on physical grounds: the real binary has more de- 
grees of freedom, reflected in the orbit and critically cur- 
rent quadrupole radiation. These differences fundamen- 
tally limit the accuracy of synthetic signals that hope to 
reproduce precessing signals using suitably-rotated non- 
precessing binaries. Equivalently, the difference between 
the red and black curves in Figure [T] suggests an accuracy 
threshold for proposals that hope to carefully calibrate 
such models against numerical merger signals. 

For comparison, we have repeated the fitting process in 
the simulation frame, computing P max by directly com- 
paring ipipB22 to ip4NR2.2- Most of our simulations 
have approximately aligned initial orbital (L) and total 
(J) angular momentum directions. As a result, the fit- 
ting procedure generally finds a similar best fit P max at a 
similar parameter location. The small offset between the 
best-fit simulation frame and corotating frame parame- 
ters can be qualitatively understood: rapid in-band pre- 
cession shifts the best-fit parameters by an amount pro- 
portional to the post-merger precession frequency. The 
question of directly comparing precessing simulations to 
nonprecessing models has considerable practical interest, 
particularly for efforts to estimate binary parameters us- 
ing simple models for the merger signal. That said, the 
results of this comparison are beyond the scope of our 
current study and not provided here. 



IV. UNDERSTANDING THE BEST-FIT 
PARAMETERS 

As described above, for each simulation and mass, 
we have found the nonprecessing model parameters 
(Msim, ?7sim, XPB,sim) such that the corotating-frame (2, 2) 
mode and the IMRPhenomB mode are most similar. As 
illustrated by Figure [2j the best-fit nonprecessing pa- 
rameters are neither constant nor trivially related to the 
precessing binary being simulated. In fact, on physi- 
cal grounds we expect and our calculations show that 
the best-fit nonprecessing model to evolve from a low- 
mass limit that reproduces qualitative features of the 
corotating-frame inspiral to a high-mass limit that re- 
produces qualitative features of the corotating-frame fi- 
nal black hole's ringdown. 



A. Low-mass limit 

In the advanced LIGO sensitive band, the gravitational 
wave signal from a low mass binary M < 3OOAf is pro- 
duced principally prior to merger, when the two objects 
can be distinguished as separate objects. At low mass 
the best-fitting nonprecessing model should resemble the 
initial binary, reproducing its orbital evolution. In other 
words, not only will the IMRPhenomB model fit, but the 
best-fitting model parameters are actually physical prop- 
erties of the binary. Specifically, the best-fitting non- 
precessing model ip4pg(XpB) should have similar mass 
ratio and "effective spin" [TU1 I3DH32"] . 

Figure [HJ compares the simulation's mass ratio 77 s ; m 
to the best-fitting IMRPhenomB mass ratio rjat for all of 
our long-duration simulations (i.e., an initial separation 
d > 7). To emphasize the inspiral, this figure adopts 
the best-fit parameters at M = 100M Q . Despite sys- 
tematic differences between the IMRPhenomB model and 
our nonprecessing simulations (green points), this figure 
and Figure [I] together suggest that both nonprecessing 
and precessing simulations and the IMRPhenomB model 
all produce nearly-indistinguishable estimates for the in- 
spiral phase of R~ 1 ip42 2 a ^ the same physical parameters. 

In practice, the short duration of most of our precess- 
ing simulations limits our ability to recover the inspiral 
signal in isolation. Equivalently, as illustrated in Figure 
[T] the best fit between IMRPhenomB and our precessing 
signals is relatively poor below 250M Q . That said, be- 
cause IMRPhenomB includes both inspiral and merger, we 
anticipate the best-fit parameters at M = 300Af Q will 
still correspond to binary parameters that approximate 
the simulated system just prior to merger, in a corotating 
frame; see the Appendix for further discussion. Figure [9] 
shows an example of the best-fitting IMRPhenomB spin 
parameter XPB,fit for a one-parameter family of simula- 
tions, evaluated at M — 300M Q . Given the fundamental 
physical differences between the IMRPhenomB model and 
simulations and given the mass used for comparison, we 
are not surprised that the best fit parameters do not lie 
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FIG. 8: Low mass limit: Recovering the mass ratio: At 

low masses and therefore early times, the best-fit corotating- 
frame model reflects the initial binary. These panels show 
comparisons between the best-fit IMRPhenomB parameter r/fit 
and the physical binary mass ratio 77, using best-fit values at 
M = 100 Mq\ compare to Figure [5] This figure explicitly ex- 
cludes all of our shorter simulations (d < 7M). Green points 
are nonprecessing spin-aligned systems; black points are long- 
duration simulations from the Sq series; and blue points are 
long-duration simulations from the Tq series. A solid blue 
line at T] a im = ^(IOOMq) is shown to guide the eye. [An 
identical color scheme is adopted in the subsequent Figure p\] 



precisely on any theoretically anticipated correlation. 

Even in the best of circumstances, we do not antici- 
pate being able to accurately reproduce all parameters 
perfectly with our fitting process. At any mass, mass 
ratio, and spin, gravitational wave signals like t/ji PB (X) 
with neighboring parameters A closely resemble one an- 
other. Depending on the particular signal model and 
parameters A, some parameter combinations are easier 
to identify than others via a fitting process. The easy-to- 
measure and hard-to-measure parameters can be iden- 
tified via a Fisher matrix [34] ; as an example, in the 
interpretation of low-mass nonprecessing binaries, a par- 
ticular correlated combination of mass, mass ratio, and 
spin is exceptionally difficult to constrain observationally 
[34) . The accuracy of our fit is also limited by by system- 
atic differences between IMRPhenomB and our signal. In 
this proof-of-concept study we will not provide a detailed 
analysis of the systematic uncertainties associated with 
this fitting procedure. 



B. High- mass limit 



FIG. 9: Low mass limit: Recovering the effective 
spin: At low masses and therefore early times, the best-fit 
corotating-frame model reflects the initial binary. These pan- 
els show comparisons between the best-fit IMRPhenomB pa- 
rameter Xeff an d Xpb [Eq. pi], using best-fit values at 
M = 300Mq. Colors indicate the Sq series (black); the Tq 
series (blue); and the T series (red). Colors indicate the Sq se- 
ries (black); precessing binaries in the Tq series (blue); aligned 
spin binaries from the Tq series (green); and the short simu- 
lations of the T series (red), Eq series (purple), and Lq series 
(dark red). A solid blue line at XPB.sim = »?PB,fit(lOOM0) 
is shown to guide the eye. Though provided at high mass 
(M = 300Mq) to insure short simulations are well-resolved, 
a similar distribution is recovered when using earlier epochs 
of longer signals (i.e., M ~ 100M©). 



nal black hole mass and spins derived from our simu- 
lations (Mf jS j m , af ]S i m ) and derived from the best-fitting 
IMRPhenomB parameters (Mf^af.fit)- 

The IMRPhenomB parameters Xsimi^sim are not trans- 
parently related to the final quasinormal mode frequen- 
cies and hence to the assumed final black hole mass and 
spin. To transform the best-fitting parameters rjat, Xpb, fit 
into a specific final black hole mass Aff fit and spin df^t) 
we adopt relationships between the "initial" and post- 
merger state that reproduce the correlations seen in our 
simulations. 

The transformation from nonprecessing binaries to fi- 
nal black holes is degenerate: the initial state is speci- 
fied by three parameters (mass, mass ratio, and at least 
one "typical" spin 10 ) while the final black hole has only 
two (mass and spin). Expressions exist which relate the 
pre- and post-merger parameters fTT], [35 39, 39J. For 
generic binaries, because spins precess between forma- 



The late-time gravitational wave signal from a merg- 
ing binary should be dominated by its leading-order 
quadrupole quasinormal mode, both in the simulation 
and corotating frame. On physical grounds, we therefore 
expect the best-fitting nonprecessing model to predict a 
similar final resonant frequency and decay timescale - 
or, equivalently, a similar final black hole mass M{ and 
spin af. As before, we will distinguish between the fi- 



10 Physically, a nonprecessing binary has four parameters: mass, 
mass ratio, and two spins. Due to symmetry, a nonprecessing 
binary's i/>4 2 2 must be reflection- and exchange-symmetric, im- 
plying it cannot trivially (e.g., linearly) depend on an antisym- 
metric combination of the binary's spins. To an excellent ap- 
proximation, the leading order quadrupole ('042 2) produced by 
a nonprecessing binary depends on only three parameters: mass, 
mass ratio, and an effective spin. 
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tion and merger, the precise relationship between the 
pre- and post-merger state depends on precisely when 
the initial binary's parameters are specified [331 SOI S3] • 
In this work, however, these approximate analytic re- 
lationships will only be applied to nonprecessing initial 
parameters. 11 In that case, the final black hole mass can 
be estimated using the following combination the pre- 
merger mass, mass ratio, and spins using [39 



E 



rad 



AI 



l-M f /M=[l-E ISCO (a)}p 



(11a) 



+ 4^ 2 [4p + 16pio(a + 1) + E lsco {a) - 1] (lib) 

(lie) 

(lid) 



Eisco = v /l -2/3r /sc . (a) 

_ _ L ■ (Si + S 2 ) 

a= m 

[p ,pi)~ (0.04827,0.01707) 



(lie) 



where rjgco ( a ) is the (dimensionless) radius of the Kerr 
horizon in Boycr-Lindquist coordinates. This approxima- 
tion assumes the final black hole mass can depend on the 
spins only through a - specifically, only through the to- 
tal spin, projected along the orbital angular momentum; 
other approximations allow for more degrees of freedom 
and physics |llj . Similarly, the final black hole dimen- 
sionless spin can be computed from the pre-merger mass, 
mass ratio, and spins using either Eqs. (6-7) from Lousto 
et al. [H] 12 or the simpler Eq. (6) from [32] : 



Of = J f /M 2 = 



\ai + d 2 q 2 + L\e\q\ 
(l + q) 2 



(12a) 



= 2a/3 + t 2 r) + t 3 ri + 



S4|Oi + a 2 q 



2|2 



(1 + g 2 ) 2 
S5V + to 2 +2 L-(d 1+q "a 2 ) (12b) 



l+q 
(*o,*3,*2) = (-2.89,2.57,-3.52) 

( S4: ,s t ) = (-0.1229,0.45) 



(12c) 
(12d) 



These expressions accurately reproduce the results of our 
nonprecessing simulations. For example, given the initial 
spins, these expressions reliably reproduce the final black 
hole mass we derive from the final horizon to significantly 
better than 1% in all well-resolved simulations. 



11 The expressions provided in the text have been extremely accu- 
rately tuned for aligned-spin binaries. Though significantly less 
accurate for precessing binaries, these expressions are a quali- 
tatively consistent relationship between our simulations' initial 
and final states when the binary has two precessing spins. In 
particular, these expressions usually predict the final black hole 
mass Mf a ; m /Af s ; m for generic precessing binaries to well within 
1% of their true value. This scale can be helpfully compared with 
the significant scatter visible in the top panel of in Figure 1101 

12 To implement this expression exactly requires knowledge of the 
spins just prior to plunge and merger. We only apply this expres- 
sion to nonprecessing systems - in fact, only to the IMRPhenomB 
best-fit parameters. 



In terms of these relationships, we derive the final 
black hole mass and spin implied by the best-fitting 
IMRPhenomB parameters -Mgt,XPB,at)?7fit by direct sub- 
stitution. For example, we evaluate the final black hole 
mass Mffix by assuming Si — XPB.fitTOi-t and S2 = 
XPBfit m 2 L, where m\^ 2 are derived from ?7fi t , in the ex- 
pression above: 



M f , fit /M fit = 1 - 77fi t (l - E lsco )(a) 



(13) 



As shown by example in Figure [lOj the best-fit pa- 
rameters Aff fit,Of,fit should be close to the final black 
hole's state. This figure compares the final black hole 
state identified in each simulation with the final state im- 
plied by the best-fitting IMRPhenomB parameters. At the 
late times and high masses studied here, both the signal 
model and IMRPhenomB will be dominated by late-time 
quasinormal mode decay, with some characteristic fre- 
quency u) 2t2 (Mf,at). Qualitatively speaking, we expect 
the fitting process enforces near-equality between uj 2i2 de- 
rived from V'4ps ancl the late time ringdown frequency 
in the corotating frame. The functional dependence of 
uj 2 2 on Mf , af determines the key features seen in Figure 
First, the mass ratio Imcj2,2/R-ew2,2 depends strongly 
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on the final black hole spin af, and, therefore, this fit- 
ting procedure should reliably determine the final black 
hole spin parameter; see, for example, the bottom panel 
in Figure |10[ Second, because of the extremely limited 
number of cycles available to constrain the final black 
hole's properties, the total black hole mass cannot be re- 
liably measured. 13 In particular, because relatively little 
mass is loss to infinity (Mf jM ~ 1), the fraction of mass 
lost to infinity cannot be reliably determined; see the top 
panel of Figure [TO] 

Above and beyond the fundamental limitations set by 
the functional form of ^2.2, the best-fit IMRPhenomB final 
spin differs systematically from the physical spin param- 
eter for three reasons. Systematic differences exist be- 
tween the IMRPhenomB model and our simulations, even 
for aligned spins; see the red curves in Figure [T] These 
differences arise due principally due to physical differ- 
ences between the IMRPhenomB model and our merger 
simulations. 14 In addition, the IMRPhenomB approxima- 
tion uses its own implicit relationship between initial 
binary and final black hole, encoded in the final ring- 
down frequency [TU]. As indicated in Figure 10 using 
red points, comparable differences between the physi- 
cal final state and a final state derived from the best-fit 
IMRPhenomB model appear when this procedure is applied 



13 A more detailed treatment of parameter estimation uncertainties 
in the ringdown signal is beyond the scope of this paper. 

14 We have simulated identical nonprecessing systems at different 
resolutions and compared them with IMRPhenomB; all results for 
these simulations are included in Figure [TJ The best-fit match 
and recovered parameters depends only weakly on resolution; by 
contrast, significant differences exist between IMRPhenomB and 
any of our models. 
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both to precessing and nonprecessing binaries. We have 
eliminated physical degrees of freedom and introduced 
systematic error by requiring the final black hole mass 
and spin depend only on the black hole spins through 
Xpb- By comparison, neither expressions we use to re- 
late the initial and final state [Eqs. (11 12)] depend on 
spins through Xpb! Instead, both expressions used here 
depend on the total spin Si + S2- Alternatively, Ta- 
ble [I] has several examples of simulations with identical 
mass ratio and %pb but producing significantly differ- 
ent af, see, for example, Sq(4,0.6,90) and Sq(4,0.6,270). 
Finally, the corotating-frame transformation involves a 
rotation at an appreciable proportion of the final black 
hole's quasinormal mode frequency [14]. As a result, we 
fully expect the best-fit parameters to be offset from the 
physical parameters of the final black hole. 



V. SYNTHETIC SIGNALS AND MISSING 
PHYSICS 

A. Synthetic signals 

As described in previous work on corotating frames 
[THHH], the similarity between corotating-frame modes 
and nonprecessing systems suggests a simple strategy 
for synthetic waveform generation. In this procedure, 
we join some nonprecessing set of modes i\)^ m to a 

time-dependent rotation R, constructively generating a 
synthetic waveform for each line of sight. For exam- 
ple, each of the corotating-frame modes could be inter- 
polated phenomenologically, starting with our collection 
of corotating-frame modes. Likewise, prior to merger 
the rotation operation could track expected trajectory 
based on adiabatic quasicircular spin-orbit evolution; af- 
ter merger, the rotation operation could be measured 
and calibrated to mimic suitable precession. This very 
aggressive strategy requires careful, complicated interpo- 
lations: each mode and spins depends on several param- 
eters (two masses and two spin vectors). Less aggres- 
sive strategies could adopt simpler models for precession 
or the inspiral signal. For example, a well-motivated, 
parameterized inspiral-merger-ringdown model like the 
effective-one-body model could track the orbital phase. 

To assess how well this strategy could perform in opti- 
mal circumstances, we assume the physical rotation and 
best-fitting IMRPhenomB can be determined for each can- 
didate physical system. In this most optimistic case, the 



match P ma x,sim [Eq. {TTj)] quantifies the degree of simi- 
larity between the simulation- frame (2, 2) mode i/j^nr 2 2 
and the (2, 2) mode generated by this synthetic proce- 
dure: [R~ l 4>i PB {\pB)]2,2- The top panel of Figure M 



shows the distribution of P, 



max,corot ■ 



In short, this strat- 



egy will be strikingly successful, fitting the IMR signal 
from our generic mergers to better than a few percent. 

More striking still, our simulations generally show 
-Pmax.corot - fmax.sim- In other words, we only needed 
to assess how similar the corotating-frame signal is to 
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FIG. 10: High mass: Best fit approximately recov- 
ers the final state: At M ~ IOOO-Mq, a scatter plot of 
the true {M B i m ,a s im) and best-fit (Mat, Oat) final black hole 
masses (bottom panel) and spins (top panel), derived from 
our precessing simulations (blue) and nonprecessing simula- 
tions (red). The best-fit final black hole properties are derived 
from the IMRPhenomB best-fit parameters using Eqs. ( 11|12 1. 
In both figures, a dotted black line is provided to guide the eye 
onto af, s i m = af,fint and similarly. Because the IMRPhenomB 
model claims to reproduce nonprecessing binaries, the dis- 
tribution of the red points can be used to estimate system- 
atic uncertainties in this method. For example, in the top 
panel, the systematic offset between the red points and the 
line a,f t aim — «f,fit suggests significant systematic differences 
between our simulations' ringdown frequencies and the model 
of IMRPhenomB. In the bottom panel, the significant scatter 
in both red and blue points reflects our inability to measure 
total masses to better than a few percent by a fitting pro- 
cedure; see also the bottom panel in Figure [2] This figure 
includes a separate point for each individual simulation, in- 
cluding some simulations with physically identical parameters 
but performed with different resolutions. To better resolve 
the final black hole quasinormal modes, we have explicitly 
excluded all simulations performed with h = M/77. 



a nonprecessing signal model, to determine how well a 
synthetic precessing search will perform. 

Serious challenges exist before this synthetic waveform 
can be implemented in practice. For example, the rela- 
tionship between the physical parameters of the binary, 
the physical rotation versus time, and the best-fitting 
nonprecessing model parameters must be carefully tab- 
ulated before these synthetic signals could be applied to 
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parameter estimation. Also, this model omits physical 
degrees of freedom and cannot reproduce the signal per- 
fectly, even with the best-possible nonprecessing signal 
model. As a practical matter, the latter limit sets the 
accuracy to which this program should be pursued in de- 
tail. 

Our study also allows us to better understand the es- 
sential features a more complete interpolation strategy 
must have. Because the best-fitting nonprecessing sys- 
tems are different at low and high mass, a robust strat- 
egy must produce waveforms that resemble different non- 
precessing systems early and late in the signal. This 
procedure could work physically, with a well-motivated 
IMR model; somewhat phenomenologically, using two 
basis epochs and interpolating the junction conditions 
versus spins; or completely phenomenologically, interpo- 
lating the whole waveform. 



B. Limitations of nonprecessing analogs and 
synthetic precessing waveforms 

Precessing binaries break a symmetry: reflection 
through the orbital plane. As a result, precessing bi- 
naries have waveforms with fundamentally more com- 
plexity, even in the corotating frame. No nonprecessing 
waveform, no matter how rotated, can reproduce that 
symmetry breaking; see Appendices |A4| and |A3| As a 
result, the straightforward approach described above has 
limited accuracy, applied to generic precessing systems; 
see the executive summary and Figure [3J In this sec- 
tion, we describe the physics missed when one assumes 
the corotating-frame waveform has an instantaneous re- 
flection symmetry. 

The straightforward approach approach described 
above has not included the direct feedback of precessing 
spins on the waveform. For example, during the inspi- 
ral the orbiting spins source current quadrupole radia- 
tion. By breaking reflection symmetry about the instan- 
taneous orbital plane, the current quadrupole and sim- 
ilar precession-sourced terms introduce behavior which 
cannot be reproduced by nonprecessing waveforms, no 
matter how rotated. As described in detail in the exec- 
utive summary, Figure [5] shows the reflection-symmetric 
and reflection-asymmetric parts of the I — 2 corotating- 
frame waveform, illustrating the small but increasingly 
important impact of symmetry-breaking terms. Prior to 
merger, the symmetry-breaking terms are small, typically 
several PN orders smaller than the dominant term for 
each mode; see Appendix | A 4| for a symmetry decompo- 
sition of the PN expansion. During and after merger, 
however, these symmetry-breaking terms become a sig- 
nificant fraction of the overall amplitude \ip4\- These 
symmetry-breaking terms reflect an instantaneous bias 
towards preferentially left- or right-handed emission |14j . 
In Figure [5] as in most cases, the dominant asymmetry 
occurs between the (2, ±2) modes. 

The effect of precessing spins and bias towards one he- 



licity or another is far from academic, particularly at high 
mass. Using a data-analysis-motivated comparison in the 
simulation frame, O'Shaughnessy et al. P3] have already 
showed that this oscillating bias towards one handedness 
or the other leads to a preferred handedness in the de- 
tected signal, changing slowly as a function of mass and 
line of sight. These modulations distort the 2, 2 mode 
and generally cannot be perfectly reproduced by a non- 
precessing waveform. Figure [5] and comparable calcu- 
lations for generic sources suggests that strong (tens of 
percent) symmetry breaking occurs ubiquitously in pre- 
cessing mergers. 

The synthetic procedure described above also neglects 
higher harmonics. Higher harmonics are well-known 
to produce observationally-accessible modulations of the 
gravitational wave signal along generic lines of sight for 
generic high-mass mergers [25] US] • As suggested by Fig- 
ure pi higher harmonics are at least as important as (and 
difficult to disentangle from) spin precession effects. Our 
calculations support previous results suggesting that the 
observationally-accessible information requires detailed 
models for several harmonics, beyond the leading order 

[SUES]. 

Our results, along with the previously mentioned helic- 
ity bias described in [14] , strongly suggest that high-mass 
precessing binaries cannot be well-understood without 
breaking reflection symmetry through the orbital plane. 
In turn, these symmetry-breaking terms cannot be un- 
derstood without modeling transverse spin dependence in 
detail. For strongly precessing asymmetric binaries, we 
do not believe the waveform generation problem cannot 
be completely decoupled from kinematics, at an observa- 
tionally relevant level. Finally, as in previous studies, we 
anticipate higher harmonics will be required to interpret 
the merger signal from generic precessing sources. 



VI. CONCLUSIONS 

In this paper we explore a simple "synthetic" model for 
the leading-order gravitational wave signal from precess- 
ing, merging binaries: a suitable nonprecessing binary, 
viewed in a suitable noninertial frame. Using a data- 
analysis-driven diagnostic, we compare the (2, 2) modes 
extracted from a large collection of binary black hole 
merger simulations to such a synthetic model, both in an 
inertial frame and a "corotating" frame which tracks the 
evolution of the binary. In all cases explored here, we find 
that the late-time inspiral and merger signal from generic 
black hole mergers can be reasonably approximated by a 
nonprecessing binary seen in a corotating frame. More- 
over, as expected on physical grounds, at early times 
the corotating-frame signal resembles emission from a 
nonprecessing binary with physically similar parameters: 
similar mass, mass ratio, and "effective spin." Our study 
restricted attention to our numerical simulations, with- 
out analytic extrapolations at early ("hybridization") or 
late times; that said, our results suggest that followup in- 
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vestigations with suitable semi-analytic hybrids can rea- 
sonably approximate generic precessing merger signals 
for all time. 

Though this approximation neglects significant physics 
and has correspondingly limited accuracy, this "syn- 
thetic" approach allows the efficient generation of compli- 
cated, multimodal signals from generic merging sources. 
Because this procedure bootstraps experience gained 
from nonprecessing sources, we strongly recommend 
more effort be devoted to modeling, approximating, or 
hybridizing the observationally-relevant higher harmon- 
ics from generic two-spin nonprecessing binaries; to mod- 
eling precession before, during, and after merger; and 
hence to generating qualitatively realistic synthetic pre- 
cessing merger signals. 

While a qualitatively adequate zeroth approximation, 
this procedure does not easily generalize to a high- 
precision quantitative approximation, with controlled er- 
ror estimates. Nonprecessing waveforms simply cannot 
self-consistently reproduce features tied to the system's 
kinematics: the orbital phase versus time; the ringdown 
mode frequencies, set by the final black hole's mass and 
spin; the ringdown mode amplitudes, which can reflect 
spin-orbit misalignment; et cetera. As a familiar ex- 
ample, in post-Newtonian calculations, time-dependent 
spin-orbit and spin-spin terms must be included in the 
orbital phase and calculated from suitable spin preces- 
sion equations. More significantly, nonprecessing wave- 
forms have a reflection symmetry and thus cannot repro- 
duce current quadrupole or similar asymmetric radiation 
modes. Our study emphasized reproducing the principal 
I = 2 emission from merging binaries; generic asymmet- 
ric precessing systems possess several strong higher har- 
monics, many of which must be included to accurately 
reproduce even a nonprecessing source |25| . Our cal- 
culations support previous results suggesting that the 
observationally-accessible information requires detailed 
models for several harmonics, beyond the leading or- 
der (SJJ I2S]- We anticipate detailed parameter estima- 
tion of high-mass (M > 1OOM ) binary mergers will 
require detailed modeling of multiple modes of generi- 
cally processing binaries, an effort in support of which 
a considerably larger sample of generic binary merg- 
ers are required. Moreover, because of the significant 
role the merger plays in comparable-mass binary black 
hole mergers with M > 2OM , we suspect that high- 
precision parameter estimation of generic precessing low- 
mass systems will also require a model for precession dur- 
ing merger that has been carefully calibrated against nu- 
merical simulations. 

To the best of our knowledge, gravitational wave de- 
tection strategies have never been tested against generic 
merger signals that self-consistently include precession 
and ringdown. Based on the quantitative similarities 
characterized in this work, we strongly recommend that 
gravitational wave data analysis strategies for high-mass 
binaries (M > 50M Q ) be tested against simple synthetic 
precessing inspiral-merger-ringdown waveforms, gener- 



ated by combining plausible rotations and nonprecessing 
binary merger signals. 



Appendix A: Mathematical methods 
1. Extracting the corotating waveform 

Particularly early in the inspiral, the gravitational 
wave signal from merging binaries can be approximated 
by the emission from instantaneously nonprecessing bi- 
naries, slowly rotated with time as the orbital plane pro- 
cesses [HI [53J 0S] • At late times, the gravitational wave 
signal will reflect perturbations of a single black hole with 
a well-identified spin axis. In both cases and in between, 
a well-chosen instantaneous or global frame can dramat- 
ically simplify the decomposition of tjj^h^t) in terms of 
spin- weighted harmonic functions ^^ lm {€). These sim- 
plifications make it easier to distinguish physically rel- 
evant from superfluous modulations; to model emission 
and generate hybrids; and to formulate tests of general 
relativity itself. 

In this paper, we adopt a preferred direction V aligned 
with the principal axes of (£( a £M/ [1^1 - The tensor 
(Ct a Cb\) is defined by the following angular integral, act- 
ing on a symmetric tensor constructed from the rota- 
tion group generators C a acting on the asymptotic Weyl 
scalar: 



(£(a£b)) 



T.i mm ' ^m'^Hm ( lm ' | A°A) | Irn) 

fdn\ip 4 \ 2 
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where in the second line we expand ip^ = 

Y.i m '*Pii m {t) Y irn 2) { e i ( l ) ) and perform the angular 
integral. The action of the rotation group generators C a 
on basis states \lm) is well-understood, allowing us to 
re-express the tensor (£( a £ b \) as [33] : 
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where Q, m = ■s/Uj + 1) — m(m + 1). In terms of these 



15 



expressions, the orientation-averaged tensor is 

7 + Re(I 2 ) Im/ 2 Re/i' 
7 -Re(J 2 ) hali 

■*■ zz 

(A2e) 
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The dominant eigendirection V of this tensor specifies 
two of the three Euler angles needed to specify a frame: 



V = (cos a sin/3, sin a sin j3, cos /3) 



(A3) 



To determine the remaining Euler angle (7), we self- 
consistently adjoin a rotation in the plane transverse to 
this orientation, to account for the gradual buildup of 
transverse phase due to precession [20j : 
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Having specified the three Euler angles that define a 



new frame, we rotate the simulation-frame Y, 
cients of ^4 to the new, time-dependent frame 
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where R(a, ft, 7) carries the z axis to V, plus a rotation 
transverse to that direction by 7. 

All simulations of the same physical system (with the 
same tetrad normalization) will agree on ij)&{t,n). The 
choice of frame at future infinity reparameterizes the 
same results. While our choice for the preferred frame 
continues to precess during and after merger, to the ex- 
tent our simulations have so far resolved, other authors' 
preferred frames may converge to a fixed frame, aligned 
with the final total angular momentum direction Jf. The 
choice of corotating frame depends on convention [44]. 
As a result, the corotating-frame waveforms we describe 
below can differ from those extracted using other conven- 
tions, with differences increasing at late times. For the 
purposes of this paper - comparison with nonprecessing 
binaries, principally of the leading-order mode - we an- 
ticipate these differences are small. 

As shown in prior work [14] and illustrated again in 
the bottom panel of Figure [IT] the preferred orienta- 
tion V and rotation 7 evolve smoothly, usually process- 
ing around the total angular momentum direction. In 
this work, we do not construct synthetic waveforms and 
therefore do not report on a functional approximation to 
the Euler angles as a function of time. That said, using a 
frame aligned with the (initial) total angular momentum 
J to define the Euler angles, the opening angle j3 is nearly 
constant; the precession angle a evolves steadily forward, 
at a precession rate set either by spin-orbit coupling (in 
the inspiral) or the final quasinormal mode frequencies 
(during ringdown) ; see Figure 8 in [14] . We therefore ex- 
pect the rotation R(t) can be easily and reliably fit as a 
function of time. 



2. Complex overlap and mass-weighted 
comparisons 

We coherently compare the (noise-free) signal expected 
along any pair of orientations with a complex inner prod- 
uct motivated by the detector's noise power spectrum 
|19[ 126) . For our purposes, numerical relativity simula- 
tions take as inputs binary black hole parameters and 
desired line of sight (denoted by A) and return the Weyl 
scalar ip^ (t) , a complex- valued function of time evaluated 
along that line of sight. For any pair of simulations and 
lines of sight, we compare tp4 and ip± by a complex over- 
lap 
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where Sh is a detector strain noise power spectrum. In 
this and subsequent expressions we used unprimed and 
primed variables to distinguish between the two wave- 
forms being compared, involving potentially distinct pa- 
rameters A' and lines of sight n'. For simplicity and to 
avoid ambiguity, in this paper, we adopt a semi-analytic 
model for the initial LIGO sensitivity [15]. As with the 
single-detector overlap, the complex overlap can be max- 
imized over the event time and polarization (t c ,ip c ) and 
by a simple Fourier transform: 
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The overlap P max is unity for identical simulations and 
lines of sight. How different must P ma x be from 1 to be 
significant? Roughly speaking, mismatch leads to de- 
tectable effects when 1 — P > l/p 2 27J, for p the sig- 
nal amplitude. Given the expected signal amplitudes for 
the first few gravitational wave events, a nonprecessing 
analog is for practical purposes indistinguishable from a 
corotating waveform with P max > 99%. 

Unlike other authors, in this paper we only employ 
waveforms extracted from the fully simulated spacetime, 
without making any attempt to hybridize that signal onto 
a post-Newtonian precursor. At the same time, we adopt 
a data-analysis-driven comparison, driven by the unique 
bandpass of a plausible detector. This comparison re- 
quires us to adopt physical timescales (oc M) and fre- 
quency scales (oc 1/M). Because the simulation has fi- 
nite duration and dynamics with a finite frequency range, 
our simulations are physically relevant only for a specific 
range of masses. The limit at the low mass end is set by 
the simulations' initial orbital frequency, which at very 
low masses (~ IOOM0) can lie in a detector's sensitive 
band. At the high mass end, the exponential decay of 
post-merger oscillations implies that our comparisons are 
contaminated by numerical noise at above ~ IOOOMq. 
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To determine these two limits unambiguously, we 
perform the following quantitative tests. To set the 
upper mass limit, we use a one-parameter family 
of simulations k similar to Sq(4, 0.6, 270, 9), per- 
formed at different resolutions. We require the overlap 
('042,2,felV'42,2,fe')/l^42,2,/cll'/'4fc be greater than 0.97 for 
all pairs of resolutions. The largest mass for which this 
bound holds is ~ 700 M@. To set the lower mass limit, 
we explored how the overlap between a signal and a trun- 
cated copy of itself changed, depending on how much of 
the early inspiral was removed; see Appendix |B 2| and 
Figure |11| To conservatively insure the signal duration 
had less than a 1% influence on the overlap, we extremely 
conservatively limited M > 4OOM for our typical short- 
duration signals (V start = 6.2M). 

Advanced detectors can nominally be sensitive to ex- 
tremely high- mass objects M > 500M Q , if frequencies 
below 40 Hz are properly calibrated. In the text, we 
have optimistically assumed all frequencies above 5 Hz 
will be calibrated, allowing ground-based interferome- 
ters to detect and measure the properties of binaries 
with total masses up to ~ 2500M Q . Realistically, how- 
ever, ultra-low-frequency detector noise and calibration 
remains a significant challenge, owing to the compli- 
cated and historically nongaussian noise in this regime 
|46j . Given the suspension and other realistic constraints, 
advanced detectors may be calibrated only for frequen- 
cies above 10 Hz, limiting gravitational wave detection to 
<15OOM . 



3. How many degrees of freedom are eliminated in 
going to a corotating frame? 

At each time step, the corotating frame expansion uses 
the extracted waveform data to reconstruct 3 Euler an- 
gles. Two of the three Euler angles, specifying the di- 
rection of V, are reconstructed from the instantaneous 
value of (jCujCi))) and do not depend on the past his- 
tory of the binary 15 The third Euler depends weakly 
on the past history of the binary, through the minimal- 
rotation condition. By contrast, each constant-/ subspace 
has 2 x (21 + 1) real degrees of freedom in the amplitudes 
and phases of its (21 + 1) modes. For example, the I = 2 
subspace alone has 10 real degrees of freedom, while the 
set of modes I < 4 has 42. By simple parameter counting, 
a corotating frame expansion cannot eliminate as many 
degrees of freedom as can exist in the waveform. 

By contrast, nonprecessing simulations exhibit many 
extremely strong symmetries between different modes. 
First, nonprecessing modes are always chiral: 



</W/) = ° 
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15 Since the orientation tensor (^Cr a Ci > )') cannot distinguish between 
izV, we additionally require V be continuous and start with V ■ 
L > 0. 



Second, nonprecessing binaries have amplitude and phase 
conjugation symmetry: ^>4* m = (— 1) "04; - m - Third, 
during the inspiral, all modes evolve in phase, with 
arg^4; m = m$ or fc. Repeatedly corroborated empirically 
[4"7] , this fact suggests that the inspiralling binary emits 
as if a rigid body. These symmetries significantly reduce 
the number of degrees of freedom needed to specify a 
nonprecessing source. To use the I = 2 subspace as an 
example and omitting the three (constant) Euler angles 
needed to describe the system, prior to merger only one 
phase and 3 amplitudes are needed to describe the sys- 
tem, while after merger, 2 phases and three amplitudes 
are needed. 

Our corotating expansion does not impose any of these 
properties. As a concrete example, consider a ficti- 
tious source producing only (2, 2) and (2, —2) modes in 
some constant frame misaligned with the global refer- 
ence frame. This source can produce each mode indepen- 
dently, with arbitrary amplitude and (chiral but other- 
wise arbitrary) phase as a function of time. Our corotat- 
ing expansion would identify the (constant) orientation 
of that frame and the (arbitrary) functional form of the 
two basis sources. 

That said, a phase-conjugate "corotating" signal 
("04/™ = ( — 1) V*4*-m = a i,m) phis an arbitrary time- 
dependent rotation has enough parameters to fit many 
but not all sources. Counting parameters, a phase- 
conjugate source would have Z+l amplitudes and I phases 
in each constant-Z subspace. Including an arbitrary ro- 
tation lets a phase-conjugate source fit 21 + 4 degrees of 
freedom. For the I — 2 subspace, however, two degrees of 
freedom remain that this form, though relatively generic, 
cannot fit: the two conjugate- antisymmetric ("odd") mo- 
ments (ipn, m = — (— 1) "04*-m = km f° r m = 1, 2). More 
generally, a phase-conjugate corotating source will never 
fit the odd moments bi^ m in each constant-/ subspace. 
Additionally, even a phase-conjugate source can conceiv- 
ably fit different orientations to each constant-l subspace. 

To conclude, nonprecessing simulations exhibit many 
strong symmetries. While we hope any "nonprecessing 
analog" will satisfy them, our corotating frame will not 
enforce them. First, each constant-/ subspace has a con- 
sistent preferred direction at each time. Second, prior 
to merger, the corotating modes have common phase 
evolution. Finally, nonprecessing systems cannot source 
current moments and must emit symmetrically when re- 
flected through the orbital plane. 



4. Symmetry expansions in the corotating frame 

To better identify and discuss small features that spin 
precession introduces into the waveform, we sometimes 
separate the corotating- frame Weyl scalar R~ l ip^ into 
different parts, reflecting symmetries. 

One way to split 1JJ4 is the usual (/ + 1 derivatives of 
the) "mass" and "current" quadrupole moments I;. m and 
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Si. m respectively, which we define as 
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where X i>m = (-l) m Z£_ m and S hm = (-l) m S£_ m have 
the usual symmetry |48| . This operation defines two pro- 
jections tpA = V>4 + ^4 that uniquely separate ip^ into 
mass and current contributions. 

Unfortunately, even the leading-order emission from 
nonspinning binaries produces both mass and current 
moments [23, 48, 49]. The standard mass and current de- 
composition does not provide a high-precision tool to dis- 
tinguish between the nonprecessing and precession-only 
contributions. 

For our phenomenological purposes, a more produc- 
tive decomposition uses classical axial and polar par- 
ity. For notational convenience, we define this split 
on a mode-by-mode basis, using conjugation symme- 
try to define CP-odd (bi >m = -{-l) l b*_ m ) and CP- 
even (a; :TO = (— 1) a* _ m ) parts of the corotating-frame 
R~ 1 ip4, as described in the text [Eq. (|3j]. As a con- 
crete example, we provide leading-order post-Newtonian 
expressions for hi m below. 

Nonprecessing binaries are even under this transforma- 
tion, radiating symmetrically above and below their or- 
bital plane. Equivalently, on a mode-by-mode basis, each 
R-handed mode m > has a corresponding L-handed 
mode m < that radiates identically but in the oppo- 
site direction. Precessing binaries break this symmetry, 
even in the corotating frame. As noted in [14], precess- 
ing binaries show a slight bias towards either R- or L- 
handed emission at any instant, with the sign of the bias 
oscillating as the orbit changes the relative orientation 
of the spins to the binary separation. This symmetry- 
breaking bias generally persists in the corotating frame 
on a mode-by-mode basis: often b^ rn ^= 0. For this rea- 
son, a conjugation-symmetry-based diagnostic provides 
a powerful tool to identify and quantify the impact of 
precessing spins in a corotating frame. 

To illustrate the expected functional form and symme- 
try properties of different multipole orders, we provide 
selected terms from the a, b symmetry decomposition of 
h^ m , as tabulated elsewhere '23\- To highlight their sym- 
metry properties, we convert their notation to explicit 
Cartesian 3-vector operations. Specifically, in place of 
the orbital phase (including tail terms) 'J, we adopt a co- 
ordinate 3-vector f for the radial separation. To convert 
between coordinate expressions and vectors, we employ 
a reference frame x, y that corotates with the binary 16 



e± 



V2 



(x ± iy) 



-V2 



r ■ e , 



(All) 
(A12) 



Our expressions follow directly from their expressions, 
substituting i = and replacing powers of exp i\t with 
\[2f ■ e* + as needed. 

For example, the leading-order I = 2 multipolcs have 
the following form, highlighting spin dependence and 
working to v A order: 



a 2 2 = spin-independent to 0(v ) 



(A13) 

(A14) 
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45 • Lv 4 + 0(w 5 ) 

(A17) 
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where to simplify comparisons with the literature we pro- 
vide the symmetry coefficients a, b derived from the com- 
plex h (i.e., -04 = dfh = d^(a h + &'')). In these expres- 
sions, we employ the same notation as |23j : 



V = 
5 = 



X± 



7711TO2 
(mi +m 2 ) 
777,1 — fn% 

mi + m 2 

(Si/m^ ± S 2 /m 2 2 )/M 



Additionally, following 
S = 



|4"4"] we introduce 

X- + 5\+ 

M \777l 7772/ 



(A19) 

(A20) 
(A21) 

(A22) 



In the notation of _23. , a would evolve with time. 



Appendix B: Simulations 
1. Simulations 

Table 1 in [H] enumerates many of the simulations 
and groups of simulations shown in this study; we 
adopt similar notation to characterize each simulation 
(T,Tq,S,Sq,. . .). To simplify the process of identify- 
ing and distinguishing between similar simulations, we 
provide each simulation with a short descriptive string. 
While the specific interpretation of the string depends on 
the simulation type, many of our simulations are denoted 
by a string of the form X(q, a, 9, d) where q is the mass ra- 
tio; a is the typical dimensionless spin; 9 is an angle; and 
d is the initial binary separation in units of M . Initial 
data was evolved with Maya, which was used in previ- 
ous binary black hole (BBH) studies [5DH55) . The grid 
structure for each run consisted of 10 levels of refinement 
provided by CARPET [SS], a mesh refinement package for 
CACTUS 160 . Each successive level's resolution decreased 
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by a factor of 2. Sixth-order spatial finite differencing was 
used with the BSSN equations implemented with Kranc 
[51]. While the initial data, grid size, resolution bound- 
aries, and peak resolution all differ between simulation, 
each simulation a member of one of the following classes: 

• S series: In this series, two equal-mass holes are 
positioned at ±3.1M, starting with S\/m\ = —ax 
and a,2 = ah(9, 0) for several choices of 9 and a. 
The outer boundaries are located at 317.4M. Each 
simulation was performed at a resolution of M/77 
on the highest refinement level. Notably, one refine- 
ment region extends between r = 20AT to r — 80M. 
As previously [55], the outermost refinement levels 
from r = 80Af to r — 317AT are too low resolution 
(dx > 3Af ) to safely reproduce fine, high-frequency 
features in the waveform. We therefore extract on 
the refinement region bounded by [20, 80] A/ (i.e., 
at r = 40, 50, 60, 75). 

• T series: In this series (also denoted the "A se- 
ries" elsewhere), the two equal-mass holes are po- 
sitioned at ±3.1M , starting with S\/m\ = az and 
S2/1TI2 = an(9, 0) for several choices of 9 and a. 
The outer boundaries are at 31 7 AT. Each simula- 
tion was performed at a resolution of M/77 on the 
highest refinement level. Refinement boundaries 
occur at r = 40, 79Af , with dx = 0.83M in that 
refinement region. As previously 62] , we therefore 
extract on r — 40, 50, 60, 75. 

• Sq series: Dimcnsionless spin vectors S^/to 2 , are 
chosen as in the S series. Each simulation was per- 
formed at a resolution of M/140 on the highest re- 
finement level. As above, we extract only in one 
refinement region, here at r = 40, . . . 90. However, 
the grid size depended on the initial starting sepa- 
ration. 

For simulations started with initial separation d = 
9Af, the outer boundaries are at 307M . Between 
r = 32Af and r = 102M , the refinement levels have 
dx = 0.9M. 

For simulations started with d = 6.2M , the outer 
boundaries are at 409 M. Between r = 25 Af and 
102M, refinement levels have dx = 0.9M. 

• Tq series: Dimcnsionless spin vectors S^/to 2 , are 
chosen as in the T series. The outer boundaries are 
at 409. 6M. Each simulation was performed at a 
resolution of A// 120 on the highest refinement level. 
Refinement boundaries occur at r = 20, 102M, be- 
tween which dx = 1AT. We extract in this refine- 
ment region along r = 40, ... 90 and extrapolate to 
infinity in that zone. 

When both spins are in the orbital plane (9 = ir/2 
or 37r/2), the initial data adopted for these se- 
ries can be comparable to evolutions performed 
elsewhere. For example, the Sq(4,0.6,90) and 



Tq(4,0.6,270) sequences are almost identical, mod- 
ulo a small shift in initial separation. 

• V series: In this series, the two equal-mass holes 
are positioned with a separation d = 6.2, starting 
with Si/m\ = ah(0, <f>) and S^/m^ = an(6, 4> + %). 
The outer boundaries are at 317M ; the highest res- 
olution is M/77; and we extract information in the 
refinement region between 20 M . . . 80 M, in which 
dx = 0.8M. 

• Lq series: Similar to the Sq series, except using a 
wider range for the smaller black hole's spin direc- 
tion. Each simulation was performed with a reso- 
lution M/140 on the highest refinement level; the 
simulation extends to r = 409. 6M; we extract only 
in one refinement region between 40Af . . . 102AT, 
where grid spacing is dx = 0.9Af . 

• Eq series: Similar to the Lq series, except starting 
at a larger separation (d = 7M) and both spins 
are in the xy plane and antiparallel, with S\/m\ = 
0.6n(7r/2, </>). Resolutions and extraction radii are 
as described above. 

• z and zq (aligned) series: Equal-mass (z) and un- 
equal mass (zq) aligned-spin binaries, with iden- 
tical specific angular momentum S/to 2 on each 
black hole. The outer boundaries are at 317AT; 
most simulations have their finest level of refine- 
ment at h — AT/103; and refinement boundaries 
occur at r = 39A/ and r = 79Af , between which 
dx = 1.24M. 

• zU (aligned, asymmetric): Unequal-mass aligned- 
spin binaries, where the two black holes do not have 
identical specific angular momentum S/m 2 . The 
outer boundaries are at r = 409. 6Af ; each simula- 
tion's finest refinement level has h = Af/160; and 
in the interval 25 . . . 102. 4Af, where we extract in- 
formation, the grid spacing is dx — 0.8A/. 

Given the significant systematic errors inherent in any 
comparison between IMRPhenomB and numerical simula- 
tions, we did not extrapolate the waveform to infinity 
or in resolution. Instead, we performed calculations on 
a fixed resolution and fixed extraction radius. We used 
multiple resolutions and extraction radii to assess our 
(relatively small) numerical error. 



2. Extrapolation, resolution, and duration tests 

Using the ip4 extracted on each of the constant- 
simulation-radius spheres listed above, we have calcu- 
lated preferred directions and corotating-frame wave- 
forms, both on constant radial slices and using ip$ ex- 
trapolated to infinity. The preferred direction agrees al- 
most exactly between these different options. Though 
the simulation-frame and corotating-frame 1JJ4 do change 
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FIG. 11: Comparing resolutions: Top panel: A plot of 
fmax,corot versus mass for the Sq(4,0.6,270,9) simulation, per- 
formed at resolutions h = M/100 (red), M/120 (blue); M/140 
(green); Af/160 (purple); and M/180 (black). Though sig- 
nificant, the differences between resolutions are nonetheless 
smaller than the typical differences between simulations; com- 
pare to Figure [T] Bottom panel: The path of the preferred 
direction derived from simulations with different resolution, 
viewed in projection in a frame aligned with the precession 
cone (i.e., with W — J; see [14jV 



slightly with extraction radius, they evolve princi- 
pally in amplitude and in common. To use a spe- 
cific quantitative measure that is directly relevant to 
our principal result, we both evaluated normalized 

overlaps (ipi 2 ,2( r )\^2,2(r')) /iV^.aMIV^CrOI between 
(2,2) modes extracted from different radial slices and 
evaluated P max using different extraction radii. For the 
selected cases examined, the choice of extraction radius 
had minimal impact (AP max < 0.001). We therefore 
limit attention to a single extraction radius (r — 75 M). 
For selected nonprecessing and generic initial config- 
urations, we have also performed simulations at suc- 
cessively higher resolutions. As concrete example, we 
have performed the Sq(4, 0.6,90,9) simulation at resolu- 
tions M/100, M/120, M/140, M/160, M/180; see Figure 



TT] Each resolution produces a slightly different dura- 
tion, with the dominant source of error being a slowly- 
accumulating phase shift, principally accumulating prior 
to merger. After recentering our simulations, we find con- 
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FIG. 12: Signal duration and the low-mass error: Top 



panel: The best-fitting match P m 



between IMRPhenomB 



and the (2, 2) mode of an unequal-mass aligned-spin binary 
[Tq(1.75,0.2,0,10)]. Dotted curve shows -P max derived using 
the whole signal; solid curve shows P max derived after first 
truncating the signal, to mimic the results of a shorter simu- 
lation starting at d — 6.2. Comparing with Figure [T] and the 
text, this illustration shows that signal duration dominates 
our error at low mass and has little impact on our results at 
high mass. Bottom, panel: At M — lOOM© (blue and black 
arrows) and M — 250M Q (purple arrows), the best-fitting 
IMRPhenomB parameters which reproduce Tq(1.75,0.2,0,10), 
truncated to different lengths. For comparison, the red point 
indicates the corresponding physical properties expected from 
the initial data. Arrows connect longer to shorter signals, in- 
dicating the effect of truncating the signal: a significant error 
in r\ and a moderate bias in \- As illustrated by the discrep- 
ancy between the best-fit points and physical parameters, dif- 
ferences between IMRPhenomB model and our simulations also 
contribute significantly to systematic error in parameter re- 
covery, particularly in mass ratio. 



vergence consistent with our differencing order, for the 
waveform epochs most significant for our results. Our de- 
fault resolutions (e.g., M/140) provide more than enough 
precision for our most striking results: the persistent pre- 
cession of V long after merger. The lower resolutions used 
on some shorter simulations are also more than adequate 
to identify principal features during merger. 

Many of our simulations commence shortly prior to 
merger, from separations as short as d ~ 6.2M. These 
short (few hundred M) simulations cannot produce a 
gravitational wave signal from an inspiral which they do 
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not simulate. To quantify the systematic errors intro- 
duced by our choice of signal duration, we have artificially 
truncated a long, high-resolution nonprecessing signal to 
a duration comparable to a typical short signal length 
from d ~ 6.2M. As illustrated in the top panel of Fig- 
ure [l2j this truncated signal matches IMRPhenomB signif- 
icantly less well than our simulation below M ~ 5OOM 
but agrees at higher masses. Since signal duration de- 
pends sensitively on spin orientation, we believe that 
most disagreement shown in Figure [l] between our simu- 
lations' corotating-frame signals and IMRPhenomB is dom- 
inated by our limited signal duration. In other words, we 
anticipate the IMRPhenomB model is an adequate repre- 
sentation of any observationally-accessible interval of the 
corotating-frame signal. 

The signal duration of our shortest simulations also 
imposed a severe limit on our ability to reconstruct the 
pre-merger binary's properties using the best-fitting low- 
mass parameters. As a concrete example, the bottom 
panel of Figure |12| compares the physical simulation prop- 
erties (red point) with the parameters derived by fitting 
IMRPhenomB to a truncated copy of a long, aligned-spin 
simulation, using M = 100M Q . Due to systematic dif- 
ferences between IMRPhenomB and our signal, the best- 
fitting and simulation parameters are significantly offset. 



The recovered mass ratio 77 at low mass depends sen- 
sitively on the simulation duration - for the case illus- 
trated here, the best-fit mass ratio changed from q ~ 1.5 
to q ~ 1. By contrast, not only was the recovered spin 
nearly independent of the simulation duration, but also 
it could be reliably estimated using a much later phase 
of the signal (i.e., with a much higher mass); see, for ex- 
ample, the purple arrows in Figure [12] and the best-fit 
spin versus mass shown in Figure [2] For this reason, in 
the text we use the recovered spin at M = 3OOM0 to 
illustrate reasonable agreement between the IMRPhenomB 
model and typical simulations. Using this choice, we 
could include many more simulations in Figure [9j 
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TABLE I: Simulations: The first column is a key, encoding the fam- 
ily, mass ratio, black hole spin magnitude(s) \Si\/M 2 and l&l/Aff and 
alignment. The next column is the simulation resolution h/M. In a 
handful of cases used as resolution tests in the text, multiple resolutions 
of the same initial data appear in the table. The next 8 columns pro- 
vide specific initial conditions: the initial separation (r start), mass ratio 
q = mi/7712, and two component spins S 2 /M 2 relative to the total initial 
mass. The column labelled T W ave provides an estimate of the duration 
of the well-resolved f/>4 2 2 mode. The column labelled xpb evaluates xpb 
[Eq. (|9|] using the specific mass ratio and spins provided in previous 
columns. Finally, the last two columns provide the final black hole mass 
and angular momentum, derived from the late-time horizon. 
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0.2 








0.2 


1900. 


0.8 


0.912 0.909 


zU(l, 0., 0.2, 11) 


160 11 


1 

















0.05 


1400. 


0.10 


0.949 0.716 


zU(l, 0., 0.4, 11) 


160 11 


1 

















0.1 


1500. 


0.20 


0.945 0.746 


zU(l, 0., 0.6, 11) 


160 11 


1 

















0.15 


1600. 


0.30 


0.942 0.775 


zU(l, 0., 0.8, 11) 


160 11 


1 

















0.2 


1700. 


0.40 


0.937 0.802 


zU(l, 0.2, 0.4, 11) 


160 11 


1 








0.05 








0.1 


1600. 


0.30 


0.942 0.775 


zU(l, 0.2, 0.6, 11) 


160 11 


1 








0.05 








0.15 


1700. 


0.40 


0.937 0.803 


zU(l, 0.2, 0.8, 11) 


160 11 


1 








0.05 








0.2 


1700. 


0.50 


0.932 0.83 



